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ABSTRACT 

We study the age and metallicity distribution function (MDF) of metal-poor stars 
in the Milky Way halo as a function of galactocentric radius by combining N-body 
simulations and semi- analytical methods. We find that the oldest stars populate the 
innermost region, while extremely metal-poor stars are more concentrated within r < 
60 kpc. The MDF of [Fe/H]< —2 stars varies only very weakly within the central 
50 kpc, while the relative contribution of [Fe/H]< —2 stars strongly increases with r, 
varying from 16% within 7 kpc < r < 20 kpc up to > 40% for r > 20 kpc. This is 
due to the faster descent of the spatial distribution (as seen from Earth) of the more 
enriched population. This implies that the outer halo < 40 kpc is the best region to 
search for very metal-poor stars. Beyond ~ 60 kpc the density of [Fe/H]< —2 stars 
is maximum within dwarf galaxies. All these features are imprinted by a combination 
of (i) the virialization epoch of the star-forming haloes, and (ii) the metal enrichment 
history of the Milky Way environment. 
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- galaxies: evolution, stellar content - 



1 INTRODUCTION 

Very metal-poor stars ([Fe/H]< —2) represent the living fos- 
sils of the first stellar generations. Their observation is cru- 
cial, as they may provide fundamental insights on both the 
properties of the first stars and on the physical mechanisms 
governing the early stages of galaxy formation, such as feed- 
back processes. An intrinsic problem that observers have to 
face is that old, very metal-poor stars, are extremely rare in 
the solar neighborhood, comprising no more that the ^0.1% 
of the stars within a few kpc of the Sun (Beers et al. 2005). 

During the past years several surveys focused on such 
elusive stellar populations, both in the Milky Way (MW) 
halo and in nearby dwarf satellites, providing an increasing 
amount of data. At the moment, the metallicity distribution 
function (MDF) of Galactic halo stars (Beers et al. 2005) 
represents one of the most important observational con- 
straints. Indeed, it consists of 2756 halo field stars observed 
within < 20 kpc of the Sun (Beers & Christlieb 2005), cov- 
ering a huge metallicity range which spans from [Fe/H]= —2 
down to [Fe/H]= —4. Despite such a large sample, the num- 
ber 300) of extremely metal-poor stars ([Fe/H]< —3) 
turns out to be still insufiicient to put solid constraints on 
the properties of the first stars (Tumlinson 2006, Salvadori, 
Schneider & Ferrara 2007), as the MDF is a rapidly increas- 



ing function of [Fe/H]. Hence, understanding where these 
stars are preferentially located is an urgent theoretical ques- 
tion. 



Carollo et al. (2007) have recently done an accurate 
kinematic study of ^ 10.000 calibration stars of the SDSS 
Data Release 5, finding that the "outer" halo, r > 15 kpc, 
includes a larger fraction of [Fe/H]< —2 stars and peaks 
at lower metallicity than the "inner" halo (r < 15 kpc). 
This evidence poses challenging questions about the physical 
origin of such segregation and the variation of the MDF with 
galactocentric radius. 



In this study we investigate the spatial distribution of 
metal-poor halo stars by combining an high-resolution N- 
body simulation for the formation of the MW (Scannapieco 
et al. 2007), with a semi-analytical model (Salvadori, Schnei- 
der & Ferrara 2007, hereafter SFS07) that follows the stel- 
lar population history and the chemical enrichment of the 
Galaxy along its hierarchical tree, successfully reproducing 
several properties of the MW and its dwarf satellites (Sal- 
vadori, Ferrara & Schneider 2008, hereafter SFS08; Salvadori 
& Ferrara 2009). 
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2 SUMMARY OF THE MODEL 

2.1 The N-body simulation 

We briefly summarize the main features of the N-body sim- 
ulation referring to Scannapieco et al. (2007), and refer- 
ences therein, for a detailed description. We simulate^ a 
MW- analog galaxy with the G CD + code (Kawata & Gibson 
2003a) using a multi-resolution technique (Kawata & Gibson 
2003b) to achieve high resolution in the regions of interest. 
The initial conditions at 2; = 56 are constructed using the 
pubhc software GRAFIC2 (Bertschinger 2001). The highest 
resolution region is a sphere with a radius 4 times the virial 
radius of the system (i.e. the MW) at z = 0, the dark matter 
(DM) particles mass and softening length are respectively 
7.8 X 10^ M© and 540 pc. The system consists on about 10^ 
particles within rvir] its virial mass and radius are respec- 
tively Mvir — 7.7x lO^^M© and Vvir — 239 kpc, roughly con- 
sistent with the observational estimates [Mvir — 1O"'^^M0, 
r^ir = 258 kpc) of the MW (Battaglia et al. 2005). The sim- 
ulation data is output every 22 Myr between 2; = 8 — 17 and 
every 110 Myr for 2; < 8. At each output a friend-of- friend 
group finder is used to identify the virialized DM haloes 
by assuming a linking parameter b = 0.15 and a threshold 
number of particles of 50. A low-resolution simulation in- 
cluding gas physics and star formation (SF) has been used 
in order to confirm that the initial conditions will lead to a 
disk formation. While gas physics is essential to reproduce 
the spatial distribution of disk stars, an N-body approach is 
suitable to investigate the halo (and bulge) population we 
are mostly interested in. 

2.2 The semi-analytical model 

We briefly describe the basic features of the model imple- 
mented in the code GAMETE (GAlaxy MErger Tree & Evo- 
lution) referring the reader to SSF07 and SFS08 for more 
details. We trace the evolution of gas and stars inside each 
halo of the hierarchy (group of DM particles) by assuming 
the following hypotheses: (a) at the highest redshift of the 
merger tree, z = 17, the gas has a primordial composition; 
(b) stars can only form in haloes of mass Mu > Ma{z) = 

3 X 1O^M0(1 + z)-^/2 {Mh > M3o{z) = 2.89 x M^z)) 
prior to (after) reionization, here assumed to be complete 
at 2; = 6; (c) in each halo the SF rate is proportional to 
the mass of cold gas; (d) according to the critical metallicity 
scenario (Schneider et al. 2002, 2006) low-mass stars with a 
Larson Initial Mass Function form when the gas metallicity 
Z > Zcr = 10-^^^ Z0; for Z < ZcT massive Pop III stars 
form with a reference mass mpopiii = 2OOM0 . 

We describe the enrichment of gas within proto- Galactic 
haloes and diffused in the MW environment, or Galactic 
Medium (GM), by including a simple description of super- 
nova (SN) feedback. Metals and gas are assumed to be in- 
stantaneously and homogeneously mixed with the gas (im- 
plications discussed in SSF07); we assume the Instanta- 
neous Recycling Approximation (IRA, Tinsley 1980). At 
each time- step the mass of gas, metals and stars within each 
halo is equally distributed among all its DM particles, and 

We adopt a ACDM cosmological model with h = 0.71, Qoh'^ = 
0.135, ^7A = 1 - Qo, Qbh'^ = 0.0224, n = 1 and erg = 0.9. 
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Figure 1. Left upper panel: comparison between the observed 
(rectangles) and simulated (shaded regions) age-metallicity dis- 
tribution of MW stars. Rectangles show the observed relation 
for different Galactic components (Freeman &: Bland-Hawthorn 
2002): the halo (blue solid rectangles); the thin and thick disk 
(black long-dashed rectangles); the bulge (violet short-dashed 
rectangles). The colored shaded areas correspond to regions that 
include, from the darkest to the lightest, the (30, 62, 90, 100)% 
of the total number of relic stars produced in the simulation, 
M^^^ ^ Ax 10^'-' M0. Top upper and lower panels: comparison be- 
tween the cumulative MDF observed in the Galactic halo within 
^ 20 kpc of the Sun (points with Poissonian error bars. Beers 
et al. 2005) and those produced in the simulation at different 
radii (violet histograms) normalized to the number of observed 
stars. The cyan histogram shows the 7 < (r/kpc) < 20 MDF for 
the inhomogeneous mixing case. For each range of radii the num- 
bers show the percentage of —2 <[Fe/H]< —1 stars (shaded area) 
with respect to the total number of [Fe/H]< —1 stars. 



such information propagated to the next integration step. 
The same procedure is applied to metals ejected into the 
GM; as a result, the chemical composition of newly virializ- 
ing haloes depends on the enrichment level of the GM out of 
which they form. Finally, the properties of long-living metal- 
poor stars hosted by each DM particle are stored to recover 
their spatial distribution at z = 0. 

As a rough estimate of the impact of the perfect 
mixing assumption above we have also explored a simple 
case of inhomogeneous metal mixing. The latter is mod- 
eled by computing the instantaneous filling factor Q = 
C^i '^7TRl{i) /3Vmw(z)) of the metal bubbles inside the crit- 
ical MW volume Vmw{z) - 30(1 + z)~^ Mpc^ and by ran- 
domly enriching a fraction F = l — exp{—Q) of GM particles. 
Note that this only provides an upper limit on F, as a clus- 
tered system, such as the stars we are considering, will have 
a smaller filling factor. A simple Sedov- Taylor blastwave so- 
lution (Sec. 3.6 of SFS07) is used to estimate the bubble 
radii Rb{i)' This case is compared with the perfectly mixed 
one in Fig. 1 (lower panel on the left). 
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Figure 2. Average density profile of —2 <[Fe/H]< —1 (top or- 
ange histogram) and [Fe/H]< —2 stars (bottom violet histogram) 
as a function of the distance from the Earth, R. Short- (long- 
) dashed curve is the /3-model (1 + {R/ Rc)^)'^^/'^ (power-law 
R~^) best-fit to the stellar distribution. 



3 RESULTS 

The model is calibrated, i.e. the SF and SN wind efficiencies 
are fixed, by simultaneously reproducing the global proper- 
ties of the MW (stellar/gas mass and metallicity) and the 
Galactic halo MDF as in SSF07. 



3.1 The age-metallicity relation 

A first test of our model results is a comparison with the 
observed stellar age-metallicity; the results are shown in the 
upper panel of Fig. 1. The simulated distribution can be 
virtually divided into three main regions defined by the pre- 
vailing formation mode of the stars contained within them. 

The first region (age > 13 Gyr, or z > 7.5) is populated 
by old stars covering almost the entire metallicity range, 
from [Fe/H]^ -0.3 to [Fe/H]= -4.2, and which correspond 
to the observed bulge and halo components (see the rect- 
angles in the panel). These stars formed in proto-galactic 
haloes associated to high (> 2a) density fluctuations, that 
virialized during the early stages of Galaxy formation at 
z > 7.5. These first stellar generations enriched the interstel- 
lar medium (ISM) of their host galaxies up to [Fe/H]> —2 (a 
process we dub as s elf- enrichment) , quenching the formation 
of additional very metal-poor stars. At the same time, met- 
als are expelled by SN feedback in the GM, thus increasing 
its metallicity above Zcr by z = 11 and allowing long- living 
metal-poor stars to form in newly virialized haloes accret- 
ing their gas from the GM. From z — 11 to z — 7.5 coeval 
formation of metal-poor and metal-rich stars occurs in differ- 
ent objects through accretion and self-enrichment processes. 
During this epoch —3 <[Fe/H]< —2 stars are produced via 
merging of self-enriched and accreted haloes. Unlike in Scan- 
napieco et al. (2007), no [Fe/H] < —4 stars form below z = 7 
because of the GM mixing approximation assumed here. 



The second region (Age < 13 Gyr, [Fe/H]< -0.3) is 
filled by stars which almost span the entire range of ages and 
metallicities along a relation on which the iron-abundance 
increases with decreasing age {halo, bulge, thick and thin 
disk). These stars formed in accreting haloes; the minimum 
[Fe/H] value of the stellar distribution at different epochs 
reflects the iron evolution of the MW environment. Stars 
located in this area only represent < 10% of the total stellar 
mass at 2; = 0. 

The third region ([Fe/H]> —0.3) is populated by iron- 
rich stars formed in self-enriched haloes, which therefore 
span all the possible ages (bulge, thin and thick disk). The 
bulk of the stellar mass resides in this region and corresponds 
to the broad peak of the SF rate (2 < ^ < 5, see for a rep- 
resentative SF history Fig. 1 of Evoli, Salvadori & Ferrara 
2008). Note that the most iron-rich stars in the simulation 
with ages < 12 Gyr, are poorer in iron than those observed. 
This systematic effect may be a consequence of neglecting 
the contribution of SNIa to gas enrichment. However be- 
cause this study is mostly concentrated on the spatial distri- 
bution of old and iron-poor Galactic halo stars this does not 
affect the main results. Even if we consider the possible exis- 
tence of a "prompt" SNIa component, with lifetime of SNIa 
of 0.1 Gyr (Mannucci, Delia Valle & Panagia, 2006), this 
is typically longer than the evolutionary time-scale of SNII 
(< 0.03 Gyr). Therefore, the iron produced by the first SNIa 
will only marginally contribute to pollute an ISM which has 
been already largely pre-enriched by several generations of 
SNII and Pop HI stars. 

In the following we will focus on the properties of old 
metal-poor [Fe/H]< —1 stars, whose features are unaffected 
by the lack of disk formation and SNIa contribution of our 
study. By excluding from our sample all the stars residing 
at distances < 1 kpc from the Galactic plane, we remove 
possible contamination by thin/thick disk stars. 

3.2 Metallicity distribution 

In the left lower panel of Fig. 1 the Galactic halo MDF 
observed by Beers et al. (2005) is compared with the sim- 
ulated one at galactocentric distances 7 < (r/kpc) < 20. 
The agreement between the model results for a radial in- 
terval representative of the observed region and the data is 
very good. In the same panel we show for comparison the 
simulated 7 < (r/kpc) < 20 MDF for the inhomogeneous 
mixing case (cyan histogram). The two mixing prescriptions 
only yield a marginal difference, mostly concentrated in the 
range —3.5 <[Fe/H]< —3 (see Sec. 4 for the discussion). 
These results encourage us to formulate explicit predictions 
for the radial dependence of the simulated MDF. To explore 
this point, we compare the same observational data with 
theoretical MDFs obtained in different radial bins and nor- 
malized to the number of observed stars. We note that the 
all-radii MDF (left upper panel) already provides a satisfac- 
tory match of the data implying that the observational sam- 
ple obtained in the above r-range provides a good proxy for 
the all-radii one. A marginal discrepancy at the low-Fe end 
of the distribution is found when comparing the observed 
data with theoretical MDFs derived for r > 20 kpc (lower 
panels), which becomes more sensible as r increases. In con- 
clusion, we find that the MDF varies only very weakly with 
radius. On the contrary, and in agreement with recent find- 
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Figure 3. Mass distribution of metal-poor (MP) -2 <[Fe/H]< -1 (left panel), very metal-poor (VMP) -3 <[Fe/H]< -2 (middle 
panel), and extremely metal-poor (EMP) —4 <[Fe/H]< —3 (right panel) stars, in the cylindrical coordinate plane (r, |C|), normalized to 
the total MW stellar mass in the simulation M*^^ ^ 4 x 10^^ Mq. 



ings by Carollo et al. (2007) and De Lucia & Helmi (2008), 
the relative contribution of —2 <[Fe/H]< —1 stars to the 
total MDF strongly depends on r, varying from 84% within 
7 kpc< r < 20 kpc, down to < 60% for r > 20 kpc. 

To better understand these features it is helpful to ana- 
lyze Fig. 3, which offers a spatial visualization of the MDFs 
in the p — C cylindrical coordinate plane, where z is the 
rotation axis and — -\- . We show the results for 
the central 100 x 100 kpc^ region. Each panel of the fig- 
ure displays a subset of relic stars in a different metallicity 
range: metal-poor (MP), —2 <[Fe/H]< —1, very metal-poor 
(VMP), -3 <[Fe/H]< -2, extremely metal-poor (EMP), 
—4 <[Fe/H]< —3. The colors show the total mass of stars 
contained in an annulus of radial width within 1 kpc nor- 
malized to the total MW stellar mass. 

The stellar distribution closely follows the dark matter 
one, i.e. it is denser towards the center and in the 10 dwarf 
galaxies found in 50 kpc < r < 100 kpc. The radial depen- 
dence of the MPs distribution (left panel) is very steep, vary- 
ing by more than 2 orders of magnitude in the inner 50 kpc; 
into the same region instead VMP/EMP stars (middle/right 
panels) are much more uniformly distributed and exhibit a 
central core. This is further illustrated in Fig. 2, which shows 
the average density profiles of MP and VMP/EMP stars 
as a function of the distance from the Earth, R. While for 
> 10 kpc the density of MP stars closely follows a power- 
law in radius, that of VMP/EMP stars is well approx- 
imated by a /3-function, [1 + (i?/i?c)^]"^^/^ with /3 = 1.4 
and Rc — 20 kpc. It follows that the relative contribution 
to the MDF of more pristine stellar generations becomes 
gradually more important at large distances (Fig. 1). Be- 
yond r ^ 60 kpc very metal-poor stars are mostly concen- 
trated within dwarf satellites, which are clearly identified 
in Fig. 3. This is in agreement with well-known evidence 
that the MDF in dwarf galaxies is shifted towards lower- 
[Fe/H] with respect to the Galactic one (Helmi et al. 2006). 
Interestingly extremely metal-poor stars are found only in 
two dwarfs and, even in these objects, they represent a sub- 
dominant stellar population (< 13%). Beyond 60 kpc the 
number of EMPs drops implying that this population is 
more condensed within such a region. 

What determines the spatial distribution of stars with 
different [Fe/H]? In addition to the underlying structure for- 
mation governed by DM, there are two key ingredients: (i) 



the virialization epoch of the star-forming haloes, which af- 
fects the final distribution of DM and hence of stars; (ii) 
the metal enrichment history of the GM, setting the ini- 
tial Fe-abundance of the ISM in newly virializing haloes. In 
Fig. 4 we show the average formation redshift of DM parti- 
cles hosting [Fe/H]< —1 stars, (z), in the p — C plane. The 
oldest stars populate the innermost region; moreover, (z) 
gradually decreases with r = vP^+C^- Beyond ^ 30 kpc 
on average (z) < 7. As for z < 8 the GM has been already 
enriched up to [Fe/H]GM ~ — 3, extremely metal-poor stars 
become more rare in such outer regions; this explain their 
spatial condensation. Very metal-poor stars, instead, extend 
up to 100 kpc as [Fe/H]GM ~ —2 when z = 5. Finally, 
as metal-poor stars predominantly form via self-enrichment 
their spatial distribution is unaffected by the GM enrich- 
ment and it is solely determined by hierarchical history of 
collapsed structures. 

A final remark concerns dwarf satellite galaxies. Beyond 
^ 30 kpc the dwarf systems can be identified as clumps of 
high (z) against the more uniform background. All satel- 
lites found in the simulation are "classical" dwarf galaxies, 
i.e. they have L > lO^L©. Only two of them, correspond- 
ing to rare > 2a fluctuations, were hosting Pop III stars 
as they virialized and began to form stars when z > 11 and 
Zgm < Zcr. The powerful explosions following the evolution 
of Pop III stars^ caused the complete blow-away of gas and 
metals; long-living stars only form at later times when more 
pre- enriched gas is collected by the dwarfs through accretion 
and merging processes. No clear imprint of their pristine for- 
mation can be found in these galaxies, which have similar 
stellar populations (([Fe/H]) ^ — 2) of "normal" dwarf satel- 
lites. Note however that these galaxies represent the most 
massive dwarfs we found with Mh = 1 — 2 x 10^ and 
M* = 0.7 - 1.2 X 10'^ M©. The remaining 80% of dwarfs cor- 
respond to < 2a fluctuations which virialize at later epochs 
z = (6 - 8) when [Fe/HjcM > -3. The lack of [Fe/H]< -3 
stars (Helmi et al. 2006) is hence naturally explained in these 
objects which have a dark matter Mh = (1 — 7) x lO^M© and 
stellar mass content M* = (0.5 — 7.5) x lO^M© consistent 
with that of the observed dwarf spheroidal galaxies. 

t Massive Pop III stars evolve as pair-instability SN. For 
^Poplll = 2OOM0 the explosion energy is 2.7 x 10^^ erg 
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Figure 4. Average formation redshift of DM particles hosting 
[Fe/H]< —1 stars in different regions of the (p, |C|) plane. 

4 DISCUSSION 

Old, [Fe/H]< —2 stars, are intrinsically rare in the Galaxy, 
representing only < 1% (i.e. < 5 x 10^ M©) of the to- 
tal stellar mass. This makes the selection of VMP stars 
one of the major challenges of stellar surveys dovoted to 
their investigation. Our study shows that: (i) the density of 
—2 <[Fe/H]< —1 stars as a function of distance from Earth 
is very steep, following a power- law, , with 7 = 2.2; 
on the contrary (ii) the density distribution of VMP/EMP 
stars exhibits a central core, closely following a /3-function, 
[(1 + (i?/i?c)^]"^^/^ with /3 = 1.4 and Rc = 20 kpc. Hence, 
though both populations are more concentrated towards the 
center, (iii) the relative contribution of [Fe/H]< —2 stars 
increases from 16% in the inner halo (at Galactocentric dis- 
tances r < 20 kpc) to > 40% in the outer halo, in good 
agreement with the observational results by Carollo et al. 
(2007,2009)^. Our findings suggest that the outer halo be- 
tween 20 kpc ^ r ^ 40 kpc is the most promising region to 
search for VMP stars, though it is obvioulsy harder to find 
more distant stars in magnitude limited surveys. 

The spatial distribution of very /extremely metal-poor 
stars is imprinted by two physical mechanisms: (i) the nat- 
ural concentration of old stellar particles towards the more 
central region of the Galaxy; (ii) the gradual enrichment, 
through SN-driven outflows, of the MW environment out of 
which haloes virialize. 

Our analysis predicts that the MDF of [Fe/H]< —2 stars 
weakly varies through the central 50 kpc. Furthermore the 
all-radii MDF, including all the stars within rvir = 239 kpc, 
satisfactory matches the data once normalized to the num- 
ber of observed stars. Two immediate consequences arise 
from these results: (i) the observed sample, which only ex- 
plored a small portion of the Galactic halo, is a good rep- 
resentation of the complete one; (ii) the results of semi- 
analytical models, which cannot incorporate the spatial in- 
formation, can nevertheless be compared with the observa- 
tions to give reasonable approximate results. 

i The low-metallicity bias in their data sample, leading to an 
underestimate of the number of [Fe/H]> —2 stars, prevents a 
rigorous comparison between simulated and observed MDFs. 



Beyond 60 kpc the density of very metal-poor stars is 
maximum inside the "classical" (L > 10^ Lq) dwarf galax- 
ies. We find that 8 out of 10 host [Fe/H]> —3 stars only, 
in agreement with observations by Helmi et al. (2006); in 
the remaining two extremely metal-poor stars represent a 
subdominant stellar population, making up only < 13% of 
the total stellar mass. Typically these galaxies virialize at 
2; 6 - 8, have a total mass = (1 - 7) x lO^M© (i.e. 
they are < 2a fluctuations of the density field), and stellar 
mass M* = (0.5 — 7.5) x 10^ M©. Interestingly, no ultra faint 
dwarf galaxies (L < 10^ Lq) are found in our simulation, 
confirming that these newly- discovered satellites are proba- 
bly left-overs of H2 cooling mini-haloes (Salvadori & Ferrara 
2009), whose physics is not included in the present study. 

We devote the final remark to our "perfect mixing" ap- 
proximation. As the porosity increases rapidly with time 
(Q ~ 1 for z ^ 7), the MDFs from the inhomogeneous mod- 
els are mostly consistent with those derived by using this 
approximation: the MDFs differences between the two mix- 
ing prescriptions (Fig. 1, left lower panel) are smaller than 
the zLla error expected from averaging over different hier- 
archical merger histories (see Fig. 6 of SFS07). On the other 
hand, the true scatter in [Fe/H] at low metallicities and large 
radii may be substantially larger than in either of these mod- 
els, as stars are strongly clustered towards the center of our 
simulation, which will reduce the true filling factor. The full 
physical modeling of metal mixing and diffusion remains one 
of the largest uncertainties in galaxy formation, and more 
work is required before one can draw definite conclusions. 
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